#! /usr/bin/sdlbrt

finp$="KFM_Player_320x200.png"
clstr=5

open "_xpm.xpm" for output as #1
file output #1,"/* XPM */"
file output #1,"static char *[]={\"4 4 16 1 \","
file output #1,"\"0 c #000000\",\"1 c #0000AA\",\"2 c #00AA00\",\"3 c #00AAAA\",\"4 c #AA0000\",\"5 c #AA00AA\",\"6 c #AA5500\",\"7 c #AAAAAA\","
file output #1,"\"8 c #555555\",\"9 c #5555FF\",\"A c #55FF55\",\"B c #55FFFF\",\"C c #FF5555\",\"D c #FF55FF\",\"E c #FFFF55\",\"F c #FFFFFF\","
file output #1,"\"0123\",\"4567\",\"89AB\",\"CDEF\"};"
close #1
shell("convert _xpm.xpm _gif.gif && rm _xpm.xpm")

xed=320:yed=200:wdwd=1
shell("rm _t*.png")
shell("convert "+finp$+" -scale "+str$(xed*8)+"x"+str$(yed*4)+"! _t2.png")
shell("convert "+finp$+" +dither -colors 16 _t4.png")
shell("convert _t4.png +dither -posterize 8 _t5.png")
shell("convert _t5.png +dither -colors 16 _t6.png")

'shell("convert _t2.png -quantize CMY -remap _t6.png -dither Riemersma _t3.png")

shell("convert _t2.png -quantize CMY -remap _gif.gif -dither Riemersma _t3.png")

setdisplay(xed*8,yed*4,32,wdwd):paper(8^8-1):ink(0):pen(0):cls
loadimage("_t3.png",1):pasteicon(0,0,1)
shell("rm _t*.png && rm _gif.gif")

dim p1[16],k1[16]

function grayscale(ct):
  cq=ct
  bq=bitwiseand(cq,255):cq=int(cq/256)
  gq=bitwiseand(cq,255):cq=int(cq/256)
  rq=bitwiseand(cq,255):cq=int(cq/256)
  return int(((rq*30)+(gq*59)+(bq*11))/100)
  end function

for y2=0 to yed-1
  setcaption(str$(int(100*(y2/yed)))+"%")
  for x2=0 to (xed*2)-1
    for y1=0 to 3:for x1=0 to 3
      xq=y1*4+x1
      p1[xq]=point(x2*4+x1,y2*4+y1)
      next:next
    for i=0 to 15
      k1[i]=grayscale(p1[i])
      next
    for j=0 to 14:for i=0 to 14
      if k1[i]>k1[i+1] then:
        tm=k1[i]:k1[i]=k1[i+1]:k1[i+1]=tm
        tm=p1[i]:p1[i]=p1[i+1]:p1[i+1]=tm
        end if
      next:next

    if clstr=0 then:
    ink(p1[ 0]):dot(x2*4+0,y2*4+0)
    ink(p1[ 6]):dot(x2*4+1,y2*4+0)    
    ink(p1[ 9]):dot(x2*4+2,y2*4+0)
    ink(p1[15]):dot(x2*4+3,y2*4+0)
    ink(p1[ 4]):dot(x2*4+0,y2*4+1)
    ink(p1[10]):dot(x2*4+1,y2*4+1)
    ink(p1[ 3]):dot(x2*4+2,y2*4+1)
    ink(p1[13]):dot(x2*4+3,y2*4+1)
    ink(p1[ 8]):dot(x2*4+0,y2*4+2)
    ink(p1[14]):dot(x2*4+1,y2*4+2)
    ink(p1[ 1]):dot(x2*4+2,y2*4+2)
    ink(p1[ 7]):dot(x2*4+3,y2*4+2)
    ink(p1[ 2]):dot(x2*4+0,y2*4+3)
    ink(p1[12]):dot(x2*4+1,y2*4+3)
    ink(p1[ 5]):dot(x2*4+2,y2*4+3)
    ink(p1[11]):dot(x2*4+3,y2*4+3)
    end if
    
    if clstr=1 then:
    ink(p1[ 0]):dot(x2*4+0,y2*4+0)
    ink(p1[ 4]):dot(x2*4+1,y2*4+0)    
    ink(p1[ 8]):dot(x2*4+2,y2*4+0)
    ink(p1[12]):dot(x2*4+3,y2*4+0)
    ink(p1[10]):dot(x2*4+0,y2*4+1)
    ink(p1[14]):dot(x2*4+1,y2*4+1)
    ink(p1[ 2]):dot(x2*4+2,y2*4+1)
    ink(p1[ 6]):dot(x2*4+3,y2*4+1)
    ink(p1[ 1]):dot(x2*4+0,y2*4+2)
    ink(p1[ 5]):dot(x2*4+1,y2*4+2)
    ink(p1[ 9]):dot(x2*4+2,y2*4+2)
    ink(p1[13]):dot(x2*4+3,y2*4+2)
    ink(p1[11]):dot(x2*4+0,y2*4+3)
    ink(p1[15]):dot(x2*4+1,y2*4+3)
    ink(p1[ 3]):dot(x2*4+2,y2*4+3)
    ink(p1[ 7]):dot(x2*4+3,y2*4+3)
    end if

    if clstr=2 then:
    ink(p1[ 0]):dot(x2*4+0,y2*4+0)
    ink(p1[10]):dot(x2*4+1,y2*4+0)    
    ink(p1[ 1]):dot(x2*4+2,y2*4+0)
    ink(p1[11]):dot(x2*4+3,y2*4+0)
    ink(p1[ 4]):dot(x2*4+0,y2*4+1)
    ink(p1[14]):dot(x2*4+1,y2*4+1)
    ink(p1[ 5]):dot(x2*4+2,y2*4+1)
    ink(p1[15]):dot(x2*4+3,y2*4+1)
    ink(p1[ 8]):dot(x2*4+0,y2*4+2)
    ink(p1[ 2]):dot(x2*4+1,y2*4+2)
    ink(p1[ 9]):dot(x2*4+2,y2*4+2)
    ink(p1[ 3]):dot(x2*4+3,y2*4+2)
    ink(p1[12]):dot(x2*4+0,y2*4+3)
    ink(p1[ 6]):dot(x2*4+1,y2*4+3)
    ink(p1[13]):dot(x2*4+2,y2*4+3)
    ink(p1[ 7]):dot(x2*4+3,y2*4+3)
    end if

    if clstr=3 then:
    ink(p1[ 0]):dot(x2*4+0,y2*4+0)
    ink(p1[12]):dot(x2*4+1,y2*4+0)    
    ink(p1[ 3]):dot(x2*4+2,y2*4+0)
    ink(p1[15]):dot(x2*4+3,y2*4+0)
    ink(p1[ 8]):dot(x2*4+0,y2*4+1)
    ink(p1[ 4]):dot(x2*4+1,y2*4+1)
    ink(p1[11]):dot(x2*4+2,y2*4+1)
    ink(p1[ 7]):dot(x2*4+3,y2*4+1)
    ink(p1[ 2]):dot(x2*4+0,y2*4+2)
    ink(p1[14]):dot(x2*4+1,y2*4+2)
    ink(p1[ 1]):dot(x2*4+2,y2*4+2)
    ink(p1[13]):dot(x2*4+3,y2*4+2)
    ink(p1[10]):dot(x2*4+0,y2*4+3)
    ink(p1[ 6]):dot(x2*4+1,y2*4+3)
    ink(p1[ 9]):dot(x2*4+2,y2*4+3)
    ink(p1[ 5]):dot(x2*4+3,y2*4+3)
    end if

    if clstr=4 then:
    ink(p1[ 0]):dot(x2*4+0,y2*4+0)
    ink(p1[ 4]):dot(x2*4+1,y2*4+0)    
    ink(p1[ 8]):dot(x2*4+2,y2*4+0)
    ink(p1[12]):dot(x2*4+3,y2*4+0)
    ink(p1[10]):dot(x2*4+0,y2*4+1)
    ink(p1[14]):dot(x2*4+1,y2*4+1)
    ink(p1[ 2]):dot(x2*4+2,y2*4+1)
    ink(p1[ 6]):dot(x2*4+3,y2*4+1)
    ink(p1[ 0]):dot(x2*4+0,y2*4+2)
    ink(p1[ 4]):dot(x2*4+1,y2*4+2)
    ink(p1[ 8]):dot(x2*4+2,y2*4+2)
    ink(p1[12]):dot(x2*4+3,y2*4+2)
    ink(p1[10]):dot(x2*4+0,y2*4+3)
    ink(p1[14]):dot(x2*4+1,y2*4+3)
    ink(p1[ 2]):dot(x2*4+2,y2*4+3)
    ink(p1[ 6]):dot(x2*4+3,y2*4+3)
    end if

    if clstr=5 then:
    ink(p1[ 0]):dot(x2*4+0,y2*4+0)
    ink(p1[ 8]):dot(x2*4+1,y2*4+0)    
    ink(p1[ 0]):dot(x2*4+2,y2*4+0)
    ink(p1[ 8]):dot(x2*4+3,y2*4+0)
    ink(p1[ 8]):dot(x2*4+0,y2*4+1)
    ink(p1[ 0]):dot(x2*4+1,y2*4+1)
    ink(p1[ 8]):dot(x2*4+2,y2*4+1)
    ink(p1[ 0]):dot(x2*4+3,y2*4+1)
    ink(p1[ 0]):dot(x2*4+0,y2*4+2)
    ink(p1[ 8]):dot(x2*4+1,y2*4+2)
    ink(p1[ 0]):dot(x2*4+2,y2*4+2)
    ink(p1[ 8]):dot(x2*4+3,y2*4+2)
    ink(p1[ 8]):dot(x2*4+0,y2*4+3)
    ink(p1[ 0]):dot(x2*4+1,y2*4+3)
    ink(p1[ 8]):dot(x2*4+2,y2*4+3)
    ink(p1[ 0]):dot(x2*4+3,y2*4+3)
    end if


    x1=x2 mod 4:y1=y2 mod 4
    u=point(x2*4+x1,y2*4+y1)
    ink(u):dot(x2,y2)
    next:next

grab(1,0,0,xed*2,yed):saveimage("_.bmp",1):shell("convert _.bmp "+finp$+"_"+str$(clstr)+".png && rm _.bmp" )
'shell("convert "+finp$+".png _.ppm && ppmtopi1 _.ppm > "+finp$+".pi1 && rm _.ppm")

for y=yed*2 to 0 step -1: for x= 0 to xed*2
  u=point(x,y/2):ink(u):dot(x,y)
  next:next

grab(1,0,0,xed*2,yed*2):saveimage("_.bmp",1):shell("convert _.bmp "+finp$+"_"+str$(clstr)+"_b.png && rm _.bmp" )
'shell("convert "+finp$+".png _.ppm && ppmtopi1 _.ppm > "+finp$+".pi1 && rm _.ppm")


'waitkey

